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ABSTRACT 

The ability to characterize exoplanets by spectroscopy of their atmospheres requires direct imaging 
techniques to isolate planet signal from the bright stellar glare. One of the limitations with the 
direct detection of exoplanets, either with ground- or space-based coronagraphs, is pointing errors 
and other low-order wavefront aberrations. The coronagraphic detection sensitivity at the diffraction 
limit therefore depends on how well low-order aberrations upstream of the focal plane mask are 
corrected. To prevent starlight leakage at the inner working angle of a phase mask coronagraph, we 
have introduced a Lyot-based low-order wavefront sensor (LLOWFS), which senses aberrations using 
the rejected starlight diffracted at the Lyot plane. In this paper, we present the implementation, testing 
and results of LLOWFS on the Subaru Coronagraphic Extreme Adaptive Optics system (SCExAO) 
at the Subaru Telescope. 

We have controlled thirty-five Zernike modes of a H-band vector vortex coronagraph in the labora¬ 
tory and ten Zernike modes on sky with an integrator control law. We demonstrated a closed-loop 
pointing residual of 0.02 mas in the laboratory and 0.15 mas on sky for data sampled using the minimal 
2-second exposure time of the science camera. We have also integrated the LLOWFS in the visible 
high-order control loop of SCExAO, which in closed-loop operation has validated the correction of 
the non-common path pointing errors between the infrared science channel and the visible wavefront 
sensing channel with pointing residual of 0.23 mas on sky. 

Subject headings: High contrast Imaging, Low-order wavefront aberrations. Extreme adaptive optics 
systems, Coronagraphy 


1. INTRODUCTION 

One of the goals of the next generation of ground- and 
space-based missions is the direct detection and spec- 
trophotometric characterization of rocky-type exoplanets 
in the habitable zone (HZ) of a parent star. The scien¬ 
tific motivation is to study the chemical compositions 
of their atmospheres to search for biosignatures. Disen¬ 
tangling rocky-type extrasolar planets from M-type and 
solar-type star at 10 parsec requires the angular resolu¬ 
tion and sensitivity of a 30-m telescope from the ground 
and 2-4 meters telescope in space respectively. However, 
resolution alone is not sufficient enough for their detec¬ 
tion in the HZ. The direct imaging of such exoplanets 
is challenged by the ability of identifying planet signal 
above the bright stellar background at small angular sep¬ 
aration, which therefore requires high contrast imaging 
(HCI) near the diffraction limit. 

Coronagraphs are used to block the starlight and sup¬ 
press the diffraction effects of the telescope, making the 
planet signal more accessible. Small inner working angle 
(IWA) coronagraphs can reach to within the first cou¬ 
ple of Airy rings of the star. However, the exploita¬ 
tion of this region relies on the ability of efficiently 
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controlling a nd calibrating the residual low-order wave- 


front errors (Guyon et al. 


2006L These aberrations oc¬ 


curring upstream of a focal plane mask (FPM) are a 
common issue for both ground- and space-based coro¬ 
nagraphs, which result in starlight leaking around the 
coronagraphic mask. The aim of this paper is to present 
the results of a unique low-order wavefront sensor appli¬ 
cable to phase mask coronagraphs (PMCs), including the 
vortex coronagraph, with which it is tested here. 

First efforts have been made to reduce the quasi-static 
pointing aberrations at Palomar well-corrected subaper¬ 
ture (WCS), on the Hale telescope, and achieved a r esid¬ 


ual of 0.02 X/D (6 mas) with a vortex coronagraph (Ser- 


abyn et al.|20I0). The current ground-based extreme 
adaptive optics (Ex AO) instruments suc h as Gemini 
Planet Imager (CPI, Macintosh et al. at the Gem¬ 

ini Observatory and bpectro-P olarimetric High-co ntrast 
Exoplanet Research (SPHERE, Beuzit et al.|20ld ) at the 
Very Large Telescope are now predictively correcting the 
dynamic low-order wavefront aberrations. 

CPI is equipped with a 7x7 low-order Shack- 
Hartmann (SH) wavefront sensor that has demonstrated 
the corrections of the non-common path aberrations 
down to < 5 nm root mean square (RMS) for spatial 
frequencies < 3 cycles/pupil under simulated turbulence. 
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By imp lementing a Linea r Quadratic Gaussian algorithm 


(LQG, Petit et al. 2009) in the AO system, they have 
demonstrated on-sky corrections of common-path vibra¬ 
tions at 60, 120 and 180 Hz to under 1 mas per axis 
for tip-tilt residuals and a reduction of focus aberration 
down to 3 nm RMS wavefront error at the 60 Hz peak 
( |Poyneer et al. 2014| . 


SPHERE’S SAXO (SPHERE AO for eXoplanet Ob¬ 
servation) uses a 40x40 visible SH wavefront sensor and 
demonstrated an on-sky residual jitter of 11 mas with an 
integrator control ler and 9 mas with an LQG algorithm 


( Petit et'ar]|2014 ). 


I'he Subaru Oo ronagraphic ExAO (SGExAO, |Jo- 
vanovic et al.|2015 ) instrument at the Subaru Telescope, 


the Exoplanetary Ci rcumstellar Environments and Disk 


Explorer (EXGEDE,|Belikov et al. 2014 Lozi et al. 2014) 


testbed at MASA Ames and the HiKh-Contrast imag- 

ing Testbed (HCIT, Kern et al. 

2013) at JPL have im- 

plemented a coronagraphic low-order wavefront sensor 

(CLOWFS, Guyon et al. 

2009a 

), which senses the re- 

jected starlight reflected 

3y the I'PM. With the use of 


a Pha se-Induced Amplitude Apodization (PIAA, [Guyon 


2003) coronagraph, residuals < 10 ^ A/H for the tip and 


tilt modes have been demonstrated in closed-loop in the 
laboratory operation. 

However, these existing solutions are not compati¬ 
ble with the non-reflective PMCs, which are the type 
of coronagraphs that diffracts the rejected starlight in 
the post-coronagraphic p upil plane. To address this is¬ 


sue, [Singh et al. (|2014a|) have introduced the concept 
Lyot-bi 


of a Lyot-based low-order wavefront sensor (LLOWFS), 
which senses aberrations using the residual starlight re¬ 
flected by the Lyot stop. Its first implementation has 
demonstrated an open loop measurement pointing ac¬ 
curacy of 10“^ \/D at 638 nm with a F our Quadrant 
Phase Mask (FQPM, ||Rouan et al. 20001 coronagraph. 
The preliminary implementation of the lLOWFS on the 
SGExAO instrument has also demonstrated an on-sky 
closed-loop pointing accuracy of ~ 7 x 10“^ \/D 


Singh et al. 2014b 

with 

WVG, n.awet et al. 

2010) 


on-sky results of an improved version of the LLOWFS 
on the SGExAO instrument. In section]^ we remind the 
reader about the principle of the LLOWFS concept and 
its integration in the SGExAO instrument. Then, sec¬ 
tion presents the results in laboratory and on-sky for 
the configuration where the aberrations sensed by the 
LLOWFS are directly corrected by the Deformable Mir¬ 
ror (DM). Finally, sectionpresents the on-sky results 
for a second configuration where the LLOWFS is inte¬ 
grated in the ExAO loop to correct for the non-common 
path and chromatic errors between the visible wavefront 
sensor of the ExAO and the imaging wavelengths. 

2. LYOT-BASED LOW-ORDER WAVEFRONT SENSOR 
2.1. Principle 

LLOWFS is a coronagraphic wavefront sensor which is 
designed to sense the pointing errors and other low-order 
wavefront aberrations at the IWA of the PMGs. The 
coronagraphic mask at the focal plane diffracts starlight 
outside the geometrical pupil in the downstream pupil 
plane. Unlike conventional coronagraphs, the diffracted 
starlight in the re-imaged pupil plane, instead of being 


simply blocked by an opaque Lyot stop, is reflected via 
a reflective Lyot stop (RLS) towards a re-imaged focal 
plane. This reflected light is collected by a detector and 
used to measure the low-order aberrations. 

LLOWFS is a linear wavefront reconstructor that relies 
on the assumption that if the post-AO wavefront resid¬ 
uals are <C 1 radian RMS then the intensity variations 
in the reflected light are a linear combination of the low- 
order aberrations occurring upstream of the focal plane 
phase mask. 

An image Ir affected by the low-order modes i of am¬ 
plitude a = (oi, 02 • • • oin) is subtracted from a reference 
image Iq and decomposed into a linear combination on 
a base of orthonormal images Si corresponding to the 
response of the sensor to the low-order modes. So the 
difference between an image at any instant and the ref¬ 
erence follows the equation 


^ ^ C^iSi 


( 1 ) 


The measurements are then used to compute the con¬ 
trol commands via an integrator control law. 

This paper focuses on the empirical approach of the 
LLOWFS only. For a detailed theoretica l description, 
the reader may refer to the publication [Singh et al.j 
([2014a[). 


2.2. SCExAO instrument with integrated LLOWFS 

SGExAO is a versatile high contrast imaging instru¬ 
ment which features an ExAO control loop using a P yra- 
mid wavefront sensor (PyWFS, Glergeon et al. 20131 that 
provides a high and stable Strehl ratio, a speckle nulling 
routine to improve the contrast on one half of the field of 
view and a LLOWFS to stabilize the starlight behind the 
coronagraphic mask. These different wavefront sensors 
are implemented on SGExAO to address the issues that 
degrade the point spread function (PSF) quality: the 
PyWFS measures the dynamical high-order wavefront 
aberrations, speckle nulling suppresses the quasi-static 
speckles and the LLOWFS measures the coronagraphic 
leaks. This publication focuses only on the LLOWFS 
and its integration with the PyWFS. More details about 
the PyWFS and the sp eckle nulling loop can be found in 
Jovanovic et al. (2015). 

'I'he SCExAO instrument is located at the Nasmyth 
platform of the Subaru Telescope. The instrument is 
sandwiched between Subaru’s 188-actuato r adaptive op- 
tics facility (AO 188, Guyon et al.[ 2014) and HiCIAO 


(Hodapp et al. 


2008), a high-contrast coronographic im- 


ager for AO offering angular/spectral/polarization dif¬ 
ferential imaging modes. Figure shows the simpli¬ 
fied version of the optical ray path on SCExAO which 
is described as follows. A0188, using the light below 
640 nm and correcting 187 modes, stabilizes the PSF 
with a typical Strehl ratio of 30% in H band. The AO 
corrected diffraction-limited F/14 beam is then fed to 
SCExAO as an input. The beam, collimated by an off- 
axis parabola (OAP), strikes SCExAO’s 2000-actuator 
DM at the pupil plane. The beam reflected from the DM 
meets the dichroic that separates the visible light (640 - 
940 nm) from the Infrared (IR) light (940 - 2500 nm). 
The visible light is reflected towards the upper bench 
via a periscope while the IR light is transmitted to the 
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Fig. 1. — Simplified optical ray path of SCExAO. The instrument is situated at the Nasmyth platform of the Subaru Telescope and feeds 
on the beam from A0188. The output of the instrument goes to the high contrast imager, HiCIAO. SCExAO has two benches: visible 
and IR. The coronagraphic masks at the focal plane are interchangeable PMCs such as VVC, FQPM and 80PM. LLOWFS is shown on 
the IR channel simply requiring a reflective Lyot stop (RLS), relay optics and a detector. The RLS presented in the figure is the Lyot stop 
designed for the VVC. 


lower bench. The visible upper bench includes a non- 
modulated PyWFS which is capable of measuring ~ 1600 
aberrated modes with a frame rate of up to 3.6 kHz at 
^ 850 nm. The lower IR bench supports the LLOWFS 
and the speckle nulling control loop working at 1.6 fj,m. 
The bench includes a variety of coronagraphs optimized 


PIAA, Shaped pupil (Kasdin et al. 

20041), VVC, FQPM 

and eight octant phase mask (8UL 

M, 

Murakami et al. 


20101. The VVC on SCExAO is a rotating half-waveplate 


structure that has a vectorial phase spiral. There is a 25- 
/rm diameter opaque metall ic spot deposited at the center 
to mask the central defect ( Mawet et al.|2009 ). We used 
this coronagraph for the results presented in this paper. 

After the dichroic, the PIAA optics mounted in a wheel 
can be moved in or out to apodize the IR beam. At the 
focal plane, all the PMCs mentioned above sit in a wheel 
that can be adjusted in the x, y and z directions via mo¬ 
torized actuators. The on-axis starlight diffracted by the 
FPMs in a downstream re-imaged pupil plane encounters 
a pupil wheel, which sits at an angle of 6° as shown in 
fig-E This pupil wheel consists of the RLSs correspond¬ 
ing to each FPMs at the focal plane. These pupil masks 
are made by lithographing a layer of chrome on a fused 
silica disk of 1.5-mm thickness. Figure E shows an ex¬ 
ample of a RLS for the VVC coronagrapm The chrome, 
corresponding to the reflective surface in this image has 
a reflectivity of only 60% in near infrared while the rest 
is being absorbed. 

The RLS at the pupil plane blocks the diffracted 
starlight rejected outside of the geometrical pupil. This 
unused masked starlight is reflected towards a Near In¬ 
frared (NIR) detector in a re-imaged focal plane for low- 
order wavefront sensing. This detector will be referred 
to as the LLOWFS camera throughout the paper. The 
nulled coronagraphic PSF is directed towards two differ¬ 
ent NIR imaging optics via a selection of beamsplitters 
that can select the spectral content and the amount of 
flux between the two optical paths. One relayed optical 
path is towards the high frame rate internal NIR imaging 
camera and another one is towards HiCIAO. 


The LLOWFS camera and the internal NIR imaging 
camera are InGaAs CMOS detectors with a resolution of 
320 X 256 pixels, a frame rate of up to 170 Hz and a read 
out noise of 140 e“. They are used for the alignment of 
the coronagraphs as well as the testing and calibration of 
the low-order control loop either with the internal cali¬ 
bration source or directly on the sky. On the other hand, 
HiCIAO uses a HAWAII 2RG detector with a resolution 
of 2048 X 2048 pixels, a frame rate ^ 3 Hz and a read 
out noise of 15 - 30 e~. HiCIAO is a facility science 
instrument we used to perform the differential imaging 
and to collect the post-coronagraphic data during the 
on-sky operations. The advantage of having both the in¬ 
ternal NIR camera and the HiCIAO is that the former 
can be used to track the high temporal frequencies in 
the atmospheric turbulence while the latter is ideal for 
tracking the slow varying spatial frequency components 
with much better sensitivity. 

The SCExAO instrument is developed with an ulti¬ 
mate goal of being rapidl y adaptable to the future ex - 
tremely large telescopes (Guyon & Martinache 2013). 
Further details of the SCExAO instrument and its fu- 
ture capabilities are beyo nd the scope of t his pa per and 
are described in detail in Jovanovic et al. (2015). 


2.3. Deformable Mirror as a wavefront corrector and a 
turbulence generator 


The DM of SCExAO cannot only be used to control 
the aberrations up to the highest spatial frequency of 
22.5 X/D but also to inject phase errors to simulate a 
dynamical turbulence for laboratory tests. The phase 
maps injected on the DM are built using a simulated 
phase screen, which follows the Kolmogorov profile. This 
phase screen can also be filtered to mimic the effects of 
the low and high spatial frequencies under pre/post-AO 
corrections. The simulated turbulence can run in the 
background independently of the corrections injected on 
the DM by the wavefront control loops. The final com¬ 
mand sent to the DM is then the sum of the injected 
turbulence and the calculated corrections. For the turbu¬ 
lence injection, we control different parameters: strength 
(amplitude in nm RMS), wind speed (m/s) and an op- 
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tional coefficient reducing the low-spatial frequencies to 
mimic the effect of the A0188 wavefront residuals. How¬ 
ever, this simulation is limited by the spatial frequency 
of the DM, which is 22.5 cycles/aperture. 

2.4. LLOWFS operation on SCExAO 

SCExAO has a dedicated low-order wavefront correc¬ 
tion loop, which uses the measurement of the LLOWFS 
to calculate the control commands. The measured aber¬ 
rations are compensated by actuating the DM by the 
following two methods: 

• Direct interaction with the DM: The low-order 
wavefront corrections are sent directly to the DM. 
In this case, 35 Zernike modes in the laboratory and 
10 Zernike modes on sky can be controlled thus far. 
The method and the results obtained are described 
in detail in section [S] 

• Indirect interaction with the DM: The second av¬ 
enue of communication is when the LLOWFS con¬ 
trols the piezo-driven tip-tilt mount of the dichroic, 
which separates visible and IR channels, to offset 
the zero-point of the PyWFS. With this configura¬ 
tion, the axis of the PyWFS is changed by mov¬ 
ing the dichroic in tip-tilt with the corresponding 
amount of measured pointing residuals. This point¬ 
ing shift in the visible channel is then compensated 
by the DM in closed-loop, hence indirectly control¬ 
ling the differential pointing errors in the IR chan¬ 
nel. We demonstrate the concept and the first on- 
sky results with this preliminary setup known as 
the differential pointing system, in section 4. 

The second approach of low-order wavefront control 
is the one that will be used in the final configuration 
of SCExAO during the scientific observations. Indeed, 
the different wavefront sensors on SCExAO use the same 
DM for the wavefront correction, therefore cannot run 
simultaneously as separate units. Nevertheless, the first 
approach is still valid for coronagraphic ExAO designs 
that have a dedicated DM for the low-order correction. 

3. LOW-ORDER CORRECTION USING DIRECT 
INTERACTION WITH THE DM 

3.1. Configuration 

Figure summarizes the configuration in a simplified 
flowchart. The starlight rejected by the coronagraph is 
reflected towards the LLOWFS camera. The reflected in¬ 
tensity at any instant is then sensed at the rate of 170 Hz. 
The low-order estimations are done by first obtaining the 
response matrix, also called calibration frames. These 
frames are acquired by applying a known amplitude of 
each Zernike mode independently to the DM. The ref¬ 
erence subtracted response of the sensor is saved as a 
response matrix. The measurements are obtained using 
the Singular Value Decomposition (SVD) method, and 
used by an integrator controller to compute the correc¬ 
tions. These corrections are then sent to the DM, which 
compensates for the low-order aberrations. 

3.2. Calibration frames acquisition 

Figure presents the response of the LLOWFS to 
probe the low-order Zernike modes. These frames are 


Low-order corrections 
sent at 170 Hz 

r-LLOWFS 

I 

F/14 output beam 1 Residual starlight 

fromA0188 i reflected by RLS 

1 ! 


Internal 

NIR camera 






/ 


1 

Simulated turbulence 
injection (internal tests) 



Fig. 2.— Flowchart of the configuration when the LLOWFS is 
directly coupled to the DM as the actuator on the IR channel of 
SCExAO. The LLOWFS camera senses the starlight reflected by 
the Lyot stop and measures the low-order aberrations. Calculated 
corrections are then sent to the DM. In this configuration, we use 
a simple integrator control law. 

acquired prior to closing the control loop. In the lab¬ 
oratory, without any simulated turbulence, we apply a 
phasemap with an amplitude of 60 nm RMS for the 35 
Zernike modes separately to the DM. The effect of these 
modes on the low-order images is subtracted from the 
reference frame to calibrate the LLOWFS response to 
the low-order modes. Figure (a)(I)-(IO) shows the re¬ 
sponse matrix obtained in the laboratory for 10 Zernike 
modes only. This figure shows a clear distinction between 
the calibration frames, indicating no confusion in the re¬ 
sponse of the LLOWFS to different low-order modes. 

In a similar manner, Fig. (b)(l)-(I0) shows the on- 
sky calibration frames obtained by applying phasemaps 
with an amplitude of 60 nm RMS for the 10 Zernike mode 
on the DM while observing the science target Epsilon 
Leonis (1.5 mas RMS of tip-tilt angle on sky). These 
calibration frames were obtained with the A0188 loop 
closed. 

The on-sky response matrix looks noisier than the one 
obtained in the laboratory. It is actually dominated by 
uncorrected phase errors, since the A0188 is the only 
loop providing wavefront correction. Even if the on-sky 
signal is not as strong as in the laboratory, the modes are 
quasi-orthogonal and still can be used to close the loop. 

3.3. Measurements 

In order to characterize the performance of a low-order 
wavefront sensor for coronagraphic purpose, it is impor¬ 
tant to understand how efficiently the pointing errors are 
measured and mitigated. We analyzed the properties like 
the linear response of the sensor, the cross coupling be¬ 
tween the low-order modes and the requirement of how 
often the calibration frames should be reacquired. 

3.3.1. Linearity 

Figure presents the linearity of the sensor to the 
tip aberration studied in case of the VVC. We applied 
phasemaps of tip aberrations with amplitudes between 
± 150 nm RMS to the DM. The impact of each phasemap 
on the low-order images was recorded. Using the re¬ 
sponse matrix acquired in Fig. ^ (a), the amount of the 
tip error as well as the residualin the other modes was 
estimated through SVD. The experiment was repeated 
20 times and the plotted data is the average of the 20 
measurements acquired. 

The linearity range of the sensor is around 
150 nm RMS (from - 50 nm to 100 nm RMS) for the tip 
mode in x. The residuals of the modes tilt in y, focus, 
oblique and right astigmatism extracted through SVD 
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(a) Laboratory 

(1)Tip (2)™ 
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Fig. 3.— Response matrix for the VVC obtained (a) in the laboratory and, (b) on sky for 10 Zernike modes. Note: These frames have 
the same brightness scale. 



Fig. 4. — Linear response of the sensor to the tip aberrations 
in the case of the VVC. The Y-axis shows the measurements esti¬ 
mated for 5 modes. The residuals of tilt, focus, oblique and right 
astigmatisms are ~ 1 nm RMS within the linear range. The blue 
dash line shows the best linear fit within the linear range (from - 
50 nm to 100 nm RMS) of the sensor. Note: The plotted data is the 
average of the aberrations estimated in a set of 20 measurements. 

are ^ 1 nm RMS within the linearity range which is a tol¬ 
erable amount of cross-coupling between the modes. The 
shift in the center of the linear range towards one direc¬ 
tion could be caused by misalignments of the beam with 
respect to the FPM, or by the 25-/im metallic dot not 
being perfectly centered with the vortex half-waveplate. 
We repeated the linearity test with the rest of the modes 
and observed a similar behavior in the range of linearity 
and the shift of the zero point. 

Therefore, the stability of the reference image on 
the low-order camera dictates how often the LLOWFS 
should reacquire calibration frames. During the acquisi¬ 
tion of the calibration, if the environmental factors, such 
as temperature variation and the flexure of the instru¬ 
ments, introduce tip-tilt errors in the reference PSF, then 
the system needs to be re-calibrated. If these drifts hap¬ 
pen prior to closed-loop operation and are out of the lin¬ 
earity range, then only the PSF need to be realigned be¬ 
hind the FPM and previously acquired calibration frames 
can be reused to close the loop. However such drifts will 
not affect the closed-loop operation as the low-order cor¬ 
rection will compensate for them. 

3.3.2. Turbulence injection in the Laboratory 

All of our experiments in the laboratory are conducted 
with simulated dynamic phase errors that were applied 


on the DM. For the turbulence simulation, we chose 
150 nm RMS as the amplitude, 10 m/s as the wind speed 
and we allowed all the low-spatial frequency components 
of the turbulence to be left uncorrected mimicking the 
case with no AO correction upstream. Figure (a) is 
the visualization of a phasemap of the simulatedturbu- 
lence applied on the DM. 



(b) Corrections appilied 
by LLOWFS 



Fig. 5. — (a) The figure shows one phase map of the dynamic 
turbulence injected into the system (on the DM) and, (b) the cor¬ 
rections computed using the LLOWFS in the laboratory. During 
closed-loop operation of the LLOWFS, the final command that is 
being sent to the DM is the sum of these phasemaps. 


3.3.3. Spectral analysis in the laboratory 

For the laboratory test presented here, the low-order 
control loop is correcting 35 Zernike modes at 170 Hz, 
the frequency of the camera. The gain of the integrator 
controller is set to 0.7. We can push the gain to high 
values because the latency of the control loop is very 
low, ^1.1 frames. 

Figure^ (b) shows the correction phasemap computed 
by the LLOWFS control loop corresponding to the tur¬ 
bulence applied in Fig.[^(a). As expected, the color map 
in both images are opposite to each other, i.e. the control 
command cancels the injected turbulence. In closed-loop 
operation, the final command applied to the DM is the 
sum of these two phasemaps. 

The frequency of the LLOWFS (170 Hz) is much higher 
than the maximum frequency resolved by the minimal 
exposure time of the science detector HiClAO (< 0.5 Hz 
for an exposure time of 2 seconds). So to have a mean¬ 
ingful evaluation of the residuals in open and closed loop, 
we will analyze them in two temporal bands: 

• 0 - 0.5 Hz : corresponds to slow varying fre¬ 
quency components temporally resolved by the sci¬ 
ence camera, i.e. the dynamical contribution of the 
turbulence in the science images of HiCIAO. 
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• 0.5 - 85 Hz : corresponds to the faster motions 
resolved by the LLOWFS but averaged by the ex¬ 
posure time of the science camera, i.e. the static 
contribution of the turbulence and the vibrations 
in the science images. 

Figure [^presents a temporal measurement of the open- 
and closed-loop residuals for 35 Zernike modes. These 
measurements (red lines) are filtered by a moving aver¬ 
age of 2 seconds to match the minimal exposure time 
of HiCIAO (black lines). In closed-loop operations, the 
stability of the residuals improved noticeably for all the 
modes. 

Figure summarizes the open- and closed-loop resid¬ 
uals for all 35 Zernike modes. We obtained a reduc¬ 
tion by a factor of 30 to 500 (median of 200) on all 
the modes for the low frequencies (< 0.5 Hz), leaving 
only sub-nanometer residuals. For the higher frequencies 
(>0.5 Hz), the factor of improvement is only between 3 
and 12 (median of 5), because it is dominated by the vi¬ 
brations that are not corrected by the controller. These 
vibrations, mostly coming from the resonance at 60 Hz of 
a Stirling cooler, are introduced by mechanical motions 
of the optical elements on the bench. In fact the vibra¬ 
tions above 10 Hz are actually amplified by the overshoot 
of the controller. The pointing residuals for open- and 
closed-loop sampled at 0.5 Hz are about 10“^ XjD RMS 
(0.8 mas) and a few 10“"^ A/H RMS (0.02 mas) respec¬ 
tively. 

The high speed of the LLOWFS helps us to analyze the 
vibrations induced either by mechanical (cryo-coolers, 
motors etc.) or environmental (telescope structure due 
to wind-shaking) factors. In order to analyze the spec¬ 
tral distribution, we study the Power Spectral Density 
(PSD) of the residuals. The PSD is calculated as the 
square modulus of the Fourier transform of the residuals. 
A Welch smoothing is performed on the PSD to reduce 
the noise. Figure [^presents the PSDs of the open- and 
closed-loop data of only the tilt mode in the laboratory. 
The improvement is about two orders of magnitude at 
0.5 Hz while high frequencies > 10 Hz are slightly am¬ 
plified. 

We have yet to identify the source of vibrations occur¬ 
ring beyond 10 Hz, which are probably due to optical 
elements vibrating inside the instrument. These oscilla¬ 
tions are for now beyond the bandwidth of the low-order 
wavefront controller and therefore amplified by its over¬ 
shoot. We are currently optimizing this control loop with 
a LQG controller to correct for the vibrations of the tele¬ 
scope and the instrument. The LQG, based on a Kalman 
filter, uses the real-time low-order telemetry to calculate 
a model of the disturbance (pointing errors, turbulence 
and vibrations) and predicts the best correction to apply. 
Further discussions of LQG implementation on SGExAO 
will be the focus of a future publication. 

3.3.4. Spectral analysis on-sky 

After having tested the LLOWFS in the laboratory 
conditions, we analyzed its performance during an on- 
sky engineering run in April 2015. The results presented 
here were taken on the science target Epsilon Leonis 
{rriff = 1.23). In this case, AOI88 closed the loop on 187 
modes providing a Strehl ratio of ^ 40% (500 nm RMS 
wavefront error) in H band. The LLOWFS then closed 


the loop on this wavefront residuals at 170 Hz with 10 
Zernike modes. Since the gain of the loop is tuned manu¬ 
ally at present, a conservative gain of 0.05 is used for this 
demonstration to ensure the stability of the closed-loop 
operation. 

Figure presents the on-sky open- and closed-loop 
residuals. Similarly to Fig.[^ the results are smoothed by 
a moving average using a window of 2 seconds to match 
the minimal exposure time of HiCIAO. The improvement 
in the closed-loop residuals is visible in the on-sky data. 
However, the residuals are more disturbed by vibrations, 
hence are noisier than those collected in the labor atory. 

The same analysis as the one explained in Sec. |3.3.3[ 
i.e. separating low frequencies below 0.5 Hz resolved by 
HiCIAO and the high frequencies above 0.5 Hz averaged 
by HiCIAO, was performed on the on-sky data and is 
presented in Fig. For low frequencies, we obtained 
a reduction by a factor of 2.5 to 4.4 (median of 3.1) for 
all the modes while for the higher frequencies, closing 
the loop corrected the residuals by a factor of 1.2 only. 
This is expected due to the small gain value of the in¬ 
tegrator controller. However, we demonstrate that the 
slow varying pointing errors are reduced down to a few 
10-3 X/D RMS (0.15 mas). 

In Fig. [m we present the on-sky PSD for the open 
and closed loop for the tilt aberration only. The profile 
of the disturbance is different from the laboratory ex¬ 
periment presented in Fig. A new vibration around 
6 Hz due to the telescope structure appeared in the on- 
sky PSD. The vibration at 60 Hz was reduced because 
the Stirling cooler causing it was removed from the in¬ 
strument. Moreover, the shape of the pointing errors is 
different from the turbulence generated in the laboratory. 
Indeed, the general slope of the PSD is smaller than a 
typical Kolmogorov distribution. The amplitude of the 
variations are sometimes larger than the linear range of 
the LLOWFS (± 170 nm RMS on the wavefront), which 
causes an underestimation of the real amplitude and a 
modification of the shape of the PSD. Due to a smaller 
gain value, the LLOWFS could not correct for the vi¬ 
brations occurring beyond 0.5 Hz but the PSD shows 
a significant improvement below that frequency. Fig. 11 
summarizes the residuals in open- and closed-loop. How¬ 
ever, because of the amplitude of the variations outside 
of the linear range, the values in the figure are probably 
underestimated. 


3.4. Processed science frames 

In this section, we present the impact of the tip-tilt and 
other low-order residuals on the frames acquired by the 
SCExAO’s internal NIR camera. Figure [^presents the 
standard deviation per pixel in a cube of 1000 science 
frames (2 ms of integration time) for open and closed- 
loop, in the laboratory (Fig. 12 (a)) and on-sky for Ep¬ 
silon Leonis (Fig. [l2| (b)), AMebaran (Fig. |l2 (c)) and 
Altair (Fig. (d))~'^These images show lower standard 
deviation for close loop images (hence darker than the 
open loop images) and a better centered beam behind 
the VVC in closed loop. However for the target Epsilon 
Leonis, the coronagraph was not centered perfectly when 
the reference frame was acquired. 

These images were obtained without the correction 
of high-order modes by the PyWFS. The LLOWFS in 
closed-loop only stabilizes the beam upstream of the 
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Fig. 6.— Residuals in open and closed loop for 35 Zernike modes obtained in the laboratory with dynamic turbulence. The red lines are 
the raw residuals while the black lines are the moving average of the residuals using a 2-second window. Figure [^quantifies the open- and 
closed-loop residuals for the measurements presented here. 



Fig. 7.— Open- and closed-loop residuals for 35 Zernike modes 
corrected under the laboratory turbulence. The correction at 
low frequencies is about two orders of magnitude, leaving sub¬ 
nanometer residuals for all the modes. 

VVC, without showing any significant contrast improve- 
ment in the absence of an Ex AO loop. Therefore, the 
on-sky contrast enhancement of the VVC cannot be eval¬ 
uated with these results. 

4. LLOWFS INTEGRATION WITH THE HIGH-ORDER 
PYRAMID WAVEFRONT SENSOR 

The final goal of the LLOWFS is to work in close in¬ 
teraction with a high-order wavefront sensor like PyWFS 
to correct for the non-common path and chromatic er¬ 
rors occurring between the imaging and wavefront sens¬ 



Fig. 8.— PSD of the open- and closed-loop for the tilt aberration 
under the laboratory turbulence. Significant improvement is visible 
in closed loop operation at frequencies < 3 Hz. The vibrations 
beyond 10 Hz are amplified due to the overshoot of the controller. 

ing channels. The control of non-common path aberra¬ 
tions is essential because the PyWFS is using the visible 
light while the coronagraph uses the NIR light. Also the 
PyWFS is not sensitive enough to low-order modes, and 
leaves a part of them uncorrected. So these uncorrected 
aberrations (static and dynamic) create unwanted stel¬ 
lar leakage around the coronagraphic mask in NIR. We 
integrated the LLOWFS with PyWFS to address these 
non-common path and chromatic errors. 

4.1. Configuration 
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Fig. 9. — On-sky open and closed loop residuals for 10 Zernike modes for the science targ et Epsilon Leonis. The red lines are the raw 
residuals whereas the black lines are the moving average with a 2-second window. Figure ^[quantifies the residuals presented here. Note: 
The open loop is the post-A0188 raw residuals and the amplitude variations of the residuals are sometimes outside of the linear range of 
LLOWFS, which cause the underestimation of their measurement. 



Fig. 10. — Open- and closed-loop residuals obtained on-sky for 10 
Zernike modes. The correction provides a significant improvement 
at low frequencies but slightly amplifies the higher frequencies. 

SCExAO’s high-order PyWFS, currently under de¬ 
velopment, is capable of controlling ^ 1600 modes at 
3.6 kHz. For the results presented in this paper, we used 
an earlier version of the PyWFS running at 1.7 kHz and 
correcting only tip-tilt. 

In this preliminary setup, PyWFS is the only system 
communicating with the DM. So instead of sending com¬ 
mands to the DM, the LLOWFS uses the differential 
pointing system to offset the zero-point of the PyWFS. 
Figure presents the flowchart of the LLOWFS inte¬ 
gration inside the high-order control loop. The blue ar¬ 
row shows the first regime described earlier in Sec. |3.1| 
(cf. Fig. [^, where LLOWFS sends commands directly 
to the DIW The red arrows describes the configuration of 
the second regime where the high-order loop corrects for 



Fig. 11. — On-sky PSD of the open- and closed-loop presented 
for the tilt aberration only. A telescope vibration around 6 Hz 
appeared during the on-sky operation. In closed-loop, an improve¬ 
ment can be noticed at frequencies < 0.5 Hz. Due to the effects 
of the non-linearities in LLOWFS response, the real amplitude of 
the residuals are underestimated, causing the slope of the PSD to 
appear smaller than the one obtained in the laboratory. 

tip-tilt aberrations in visible and the low-order loop send 
commands to the differential pointing system to compen¬ 
sate for chromatic errors in IR. 

4.2. On-sky demonstration 

We pointed the telescope to the variable star x Cyg 
(rriH = —1.1 during this observation). A0188 closed the 
loop on 187 modes with a seeing of 0.8” at 1.6 /im. The 
PyWFS closed its loop only on tip-tilt in the visible with 
a 1.7 kHz loop speed. The PyWFS was not optimized at 
this point and hence provided only a partial correction 
of tip-tilt modes. 
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(a) Laboratory 



Open Loop Closed Loop 


(b) Epsilon Leonis 




Open Loop 


Closed Loop 


(c) Aldebaran 



“-OOP Closed Loop 


(d) Altair 



Open Loop Closed Loop 


Fig. 12. — Comparison of the standard deviation of the intensity for 1000 frames of the NIR camera (a) laboratory, (b) Science target 
Epsilon Leonis {mu = 1.23), c) Science target Aldebaran {mu = —2.78) and the (d) science target Altair {mu = 0.10). Note: Each set of 
open- and closed-loop images are of same brightness scale. Closed-loop images are expected to be darker than the open loop images. Black 
spot at the middle of all the frames is the metallic dot at the center of the VVC to mask its central defects. 



Fig. 13.— Flowchart of the LLOWFS functioning in two configurations on SCExAO. The black arrows depicts the common flow of the 
low-order control loop in both regimes. Confi gura tion 1 (blue arrow) is when the LLOWFS is used directly with the DM to correct for the 
low-order aberrations as presented in section [3.1| Configuration 2 (red arrows) is when LLOWFS, after sensing differential tip-tilt errors 
in IR channel, updates the zero-point of the PyWFS using a differential pointing system. To compensate for the beam shift in the visible 
channel, the high-order loop commands the DM to correct for the chromatic errors. 
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Similar to the first configuration explained in Sec. |3.H 
the LLOWFS first acquired a response matrix in order 
to measure the non-common path errors. The on-sky ref¬ 
erence is moved in x and y with an angle of 1.5 mas to 
obtain the calibration frames for differential tip and tilt. 
Using this response matrix, we closed the LLOWFS loop 
with a gain value of 0.03. A small gain was used because 
of the slow response of the piezos driver, controlled only 
up to 5 Hz. Figure [T4| shows the successful loop closure of 
the PyWFS and the LLOWFS. Once again, the data pre¬ 
sented here are smoothed to simulate an exposure time 
of 2 seconds. When PyWFS loop is closed, we see a slight 
improvement in the stability, but a significant amount of 
non-common path residuals are still visible. These differ¬ 
ential errors are improved when low-order loop is closed. 

Table summarizes the open and closed-loop residuals 
for high and low-order control loops. Once again, we ana¬ 
lyzed the data at two different spectral bands. This table 
shows that we have achieved a factor 3 to 4 improvement 
in correcting differential tip-tilt residuals with the gain 
of 0.03 for the slow varying frequencies. As expected, the 
improvement for the higher frequencies is not significant 
due to the small gain used. 

Due to the fact that the variations are larger than the 
linear range of the LLOWFS, the residuals in tableare 
probably underestimated again. Even in such circum¬ 
stances, closed-loop pointing residuals are only about 
6 X 10“^ \/D (0.23 mas) when the dataset is sampled at 
the frame rate of the science camera (0.5 Hz). 


300 

200 

100 

0 

? 

C 

nj 
3 

■g 
(/) 

0) 

QC 200 
100 


0 50 100 150 200 250 300 


Time (s) 

Fig. 14. — On-sky open- and closed-loop residuals of low-order 
control integrated in the high-order corrections of post-AO 188 
wavefront residuals. The black data is the moving average of resid¬ 
uals with 2 second window while the red data are the raw resid¬ 
uals. When the low-order loop is open then the high-order loop 
is correcting the pointing errors only in the visible leaving chro¬ 
matic errors uncorrected. These chromatic errors are significantly 
reduced when the loop is also closed using the LLOWFS. Table ^ 
summarizes low-order residuals for the differential tip-tilt modes. 
(Science target: x Cyg.) 


PyWFS closed loop+ i PyWFS closed loop+ 

LLOWFS open loop * LLOWFS closed loop 


I - Tipi 



1 


residual turbulence for low frequencies (< 0.5 Hz). In 
closed loop, an overshoot between 1 and 2 Hz is also vis¬ 
ible. This is due to the mismatch between the frequency 
of the sensor (170 Hz) and the frequency of the actu¬ 
ator (5 Hz). A gain of 0.03 actually corresponds to a 
gain of 1 at the speed of the actuator, which explains the 
overshoot. 



Frequency (Hz) 


Fig. 15. — On-sky open- and closed-loop PSD of the differential 
tip aberration in case of the PyWFS integration with the LLOWFS. 
Closing the loop wit h th e PyWFS reduces the telescope vibrations 
at 6 Hz shown in Fig. |ll| The low-order correction provides signifi¬ 
cant improvement at frequencies <0.5 Hz and an overshoot around 
1 Hz because of the difference in the sensing (170 Hz) and the 
correction (5 Hz) frequency. 



4.3. Limitations with the initial setup 

The performance of the LLOWFS with its preliminary 
integration with the PyWFS was constrained due to sev¬ 
eral factors. 

• Due to the slow response of the piezo driver (every 
0.2 seconds), the LLOWFS could not control tip- 
tilt aberrations faster than 1 Hz. In the current 
configuration, we updated the differential pointing 
system by replacing the control of the tip-tilt from 
the piezo-driven dichroic to a tip-tilt mirror which 
is used for the modulation of the PyWFS. This will 
increase the loop rate up to 100 Hz. 

• Using either the dichroic or the tip-tilt mirror, the 
low-order control is limited to only tip and tilt 
modes. To correct other low-order aberrations as 
well and to improve the speed, we are currently 
upgrading the way the LLOWFS interacts with the 
PyWFS. The LLOWFS will send its corrections di¬ 
rectly to the PyWFS that will then overwrite its 
reference point to compensate for these corrections 
with the DM. 


We present the on-sky PSD of the high and low-order 
integrated control loops for the differential tip aberration 
in Fig. Compared to the Fig. the fast high-order 
control loop has diminished the telescope vibrations pre¬ 
viously noticed at 6 Hz. When we close the loop using 
the LLOWFS, we observe a significant reduction of the 


• The PyWFS, in its initial stage corrected only tip- 
tilt in the visible, which could not provide signih- 
cant improvement for LLOWFS in the IR channel. 
Hence, the LLOWFS performance in both config¬ 
urations was dominated by the uncorrected higher 
order modes. 























Low Preq. (< 0.5 Hz) High Freq. (> 0.5 Hz) 
(Resolved in HiCIAO) (Averaged in HiCIAO) 
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Mode 

Unit 

Open loop 

Closed loop 

Open loop 

Closed loop 


nm 

26.1 

9.4 

144 

142 

Tip 

A/D 

1.6 X 10-2 

5.9 X 10-2 

9.0 X 10-2 

8.8 X 10-2 


mas 

0.66 

0.24 

3.6 

3.6 


nm 

36.3 

9.3 

170 

166 

Tilt 

A/D 

2.3 X 10-2 

5.8 X 10-2 

10.6 X 10-2 

10.4 X 10-2 


mas 

0.91 

0.23 

4.3 

4.2 


TABLE 1 


On-sky open- and closed-loop residuals of differential tip-tilt with the low-order loop integrated with the high-order loop. The correction 

is only significant for low frequencies. 


5. LLOWFS COMPATIBILITY WITH CORONAGRAPHS 


The LLOWFS is compatible with a family of small 
IWA PMCs. Similar closed-loop laboratory performance 
has been obtained for the FQPM and the 80PM coro- 
nagr aphs as demonstrated for the VVC in the sec¬ 
tion [333 Detailed closed-loop performance analysis of 
the LLOWFS with different PMCs is intended for the 
future publication. 

However, for coronagraphs such as the PIAA using 
an amplitude mask, the LLOWFS sensing capability de¬ 
pends on the size of the FPM. An amplitude mask big¬ 
ger than the PSF core blocks most of the starlight and 
diffracts only a small fraction of it in the re-imaged pupil 
plane. LLOWFS, in that case, does not get enough 
starlight photons, hence cannot provide an optimal solu¬ 
tion. However, we have closed the loop with PIAA and 
shaped pupil with an opaque binary FPM about half of 
the size of the PSF in the laboratory. For this type of 
coronagraphs, the CLOWFS would be a more efficient 
wavefront sensor. However, it requires some hardware 
changes on SCExAO that will not be compatible with 
the PMCs. 

A way to make LLOWFS also efficient with the ampli¬ 
tude masks is to use a conic-shaped FPM that diffracts 
the starlight in a ring around the pupil in the Lyot plane. 
Such a mask provides an optimal number of photons for 
LLOWFS independently of the size of the mask. We have 
tested this solution with an achromatic phase-s hifting fo¬ 
cal plane mask (AFPM, Newman et al. 20141, which is 
based on a diffractive optical filtering technique scaling 
the size of the FPM linearly with the wavelength. This 
mask has a cone structure at its center with an angle op¬ 
timized for the residual starlight to fall within the reflec¬ 
tive zone of the RLS. The testing of AFPMs with PIAA 
and shaped pupil are currently ongoing on our instru¬ 
ment and the performance of the low-order correction in 
the laboratory and on the sky will be discussed in future 
publications. 


6. CONCLUSION 

Small IWA phase mask coronagraphs, which enables 
high contrast imaging at small angular separations, are 
extremely sensitive to tip-tilt errors. It is crucial to 
decrease these effects using all the rejected starlight 
available, which is typically discarded in a corona- 
graph. Hence, to overcome the consequences of wave- 
front aberrations at/near the diffraction limit, imple¬ 
menting LLOWFS-like technology is crucial to control 
starlight leakage around the coronagraphic mask. We 


have demonstrated the first successful on-sky closed-loop 
test of low-order corrections using LLOWFS with the 
vector vortex coronagraph on the SCExAO instrument. 

Both in the laboratory and on-sky, we showed an 
improvement of the low-order slow varying residuals 
(< 0.5 Hz), dynamically resolved by the exposure time 
of the science camera HiCIAO. In the laboratory, we ob¬ 
tained a correction of about 2 orders of magnitude for 35 
Zernike modes, while on sky, due to the use of a small 
conservative gain for the controller, the improvement is 
only a factor of 3 for 10 Zernike modes. 

We also demonstrated the capacity of the low-order 
control loop to be combined with the higher-order loop 
for the correction of the non-common path and chro¬ 
matic aberrations between this high-order loop and the 
coronagraph. We obtained a factor 3 to 4 improvement 
in a preliminary setup, using a slow differential pointing 
system. These results are expected to improve with the 
better integration of the low-order differential control in 
the high-order loop. 

Corrections of higher-order modes other than just tip- 
tilt by PyWFS should provide a Strehl ratio > 90%. 
Moreover, the implementation of a LQG control law in 
the low-order correction should significantly reduce the 
coronagraphic leakage in the IR channel. Eurther per¬ 
formance testing of the integrated control loop on-sky is 
scheduled for the upcoming observational nights at the 
Subaru Telescope. A significant enhancement in the de¬ 
tection sensitivity of the SCExAO instrument is expected 
during the future science observations. 


7. FUTURE DEVELOPMENT 

Euture work related to the LLOWFS is focused on 
three areas which are envisioned to provide high contrast 
at small angular separation. The goals are: optimal con¬ 
trol of the low-order aberrations, point spread function 
cali bration close to/ne ar the IWA using low-order teleme- 


20111 and interaction between speckle 


try (Vogt et al. 

calibration and low-order control. 

We are currently in the process of implementing a LQG 
controller for the LLOWFS on SCExAO. In order to im¬ 
prove the post processing of the science images, we will 
use the low-order telemetry of the residuals left uncor¬ 
rected by the control loop to calibrate the amount of 
starlight leakage at small angular separations. We are 
also currently studying the interaction of speckle cali¬ 
bration with the LLOWES especially for the correction 
of speckles at small IWA. 

The development and the implementation of the above 
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mentioned technologies on SCExAO should significantly 
improve the contrast around the first couple of Airy disks 
of the star. Such advancements will allow SCExAO to 
detect young Jupiters (a few Mj ) by a factor of ~ 3 closer 
to their host stars than is currently possible with other 
ground-based ExAO systems. 

Our goal is to demonstrate innovative wavefront con¬ 
trol approaches that are central to future high contrast 
systems. To maximize the performance of the corona- 
graphs by efhciently controlling and calibrating the wave- 
front at the small angular separations, we aim to search 
the best instrumental parameter space to combine the 
optimized LLOWFS control with the PSF calibration 
and speckle nulling. A precursor of these approaches im¬ 
plemented on the next generation extremely large tele¬ 


scopes and future larger space missions should enable 
direct imaging and low resolution spectroscopy of Earth¬ 
like planets in the HZ of M-type and F,G,K-type nearby 
stars respectively. 
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